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ASYMPTOTIC  BEHAVIOR  OF  AN  I NTEGRO- DIFFERENTIAL  EQUATION 


by 

Jack  K.  Hale 


Abstract :  Consider  the  functional  differential  equation 

r° 

x(t)  =  -  a(-e)g(x(t+6))d0 

*  -1 

where  a  >  0,  g  are  continuous,  a(l)  =0.  If  a  >  0  and 
fx 

G(x)  -  g  00  as  |  x  |  ->  then  it  is  known  that  every  solution 
J0 

either  approaches  a  zero  of  g  as  t  -*•  00  or  has  an  w-limit  set 

which  is  a  1-periodic  solution  of  the  equation  x  +  a(0)g(x)  =  0. 

If  there  are  only  a  finite  number  of  equilibrium  points  and  1-periodic 

orbits,  then  there  is  a  maximal  compact  invariant  set  A  in 

a,g 

C([-1,0]JR)  which  is  uniformly  asymptotically  stable.  When  a  is 

convex,  the  topologically  structure  of  A  and  the  flow  on  A 

a,g  a,g 

are  discussed  as  a  function  of  g.  When  g  is  fixed,  xg(x)  >  0, 
x  /  0,  the  complete  bifurcation  diagram  is  given  for  a  in  a  neighbor¬ 
hood  of  a  linear  function. 
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1.  Introduction.  Consider  the  functional  differential  equation 

f0 

(1.1)  x(.t)  =  -  j  a(-0)g(x(t+9))dO 

where  a  >  O,  g  are  continuous,  a(l)  =0.  If  a  >  0  and  G(x)  = 
fx 

g  -►  »  as  |  x  |  ■*  ®,  then  every  solution  approaches  a  zero  of  g 

Jo 

as  t  -*•  «>.  If  a(s)  is  linear,  then  every  solution  either 

approaches  a  zero  of  g  as  t  -*■  00  or  has  an  w-limit  set  which 

is  a  1-periodic  solution  of  the  euqation  x  +  a(0)g(x)  =  0  (see 

Levin  and  Nohel  [  6  ]  ,  Hale  (  3 ]). 

With  some  other  minor  restrictions,  there  is  a  maximal  compact 

invariant  set  A„  of  Equation  (1.1)  in  the  space  C  =  C([-1,0]^R) 

a » g 

which  is  uniformly  asymptotically  stable.  The  primary  purpose  of  this 
paper  is  to  study  the  manner  in  which  A  and  the  flow  on  A 
induced  by  Equation  (1.1)  depend  upon  the  functions  a,g.  With  the 
kernel  a  fixed  convex  function,  it  appears  that  the  flow  on  A 

a  >  a 

depends  on  more  than  just  the  zeros  of  the  function  g.  For  a 
fixed  function  g,  there  are  linear  kernels  ap  where  the  topological 
structure  of  A  can  change  for  a  in  small  neighborhood  in 

S  9  a 

C([0,1]JR)  of  aQ;  that  is,  aQ  is  a  bifurcation  point.  These 
bifurcation  points  are  analyzed  in  detail  and  are  shown  to  be 
analogous  to  bifurcation  form  a  focus.  However,  regardless  of  the 
nature  of  g,  the  degeneracy  at  the  focus  is  not  of  order  one;  that 
is,  it  is  not  the  usual  generic  Hopf  bifurcation. 
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2.  Definitions  and  background  material.  If  C  =  C([-I,0],1R)  and 
x  is  a  continuous  function  on  [u-l.o+o],  a  >  0,  we  let  x^.  £  C 
be  defined  by  xt(e)  =  x(t+0),  -1  <  6  <  0,  t  £  [o,o+a).  If 
g  €  C  (3R ^IR )  ,  k  >  1,  and  a(s)  >0  is  continuous  on  [0,1],  let 
x (4> )  be  the  solution  of  the  functional  differential  equation 

f0 

(2.1)  x(t)  =  a(-0)g(x(t+0))d6 

} -1 

with  Xg  (‘f )  =  <}> .  I  f  wc  assume  that  x(<J>)(t)  for  every  <P  €  C  is 

defined  for  t  >  -1  and  let  T  (t)<f>  =  x.(<J>),  t  >  0,  then 

a ,  g  t 

T  (t)  :  C  ->  C,  is  a  strongly  continuous  semigroup  with  T  ( t ) 4> 

a  >  g  a  >  g 

having  continuous  derivatives  with  respect  to  <P  up  through  order  k. 

Definition  2.1.  For  Equation  (2.1)  a  set  B  in  C  attracts  a  set  H 
in  C  if,  for  any  e  >  0,  there  is  a  tQ  =  tQ  (H,e)  such  that 

Ta  g(t)H  is  in  an  e-neighborhood  of  B  for  t  >  tg.  For 
Equation  (2.1),  a  set  B  in  C  attracts  points  if  B  attracts 

the  set  {<t>}  for  every  4>  €  C .  Equation  (2.1)  is  point  dissipative 

if  there  is  a  bounded  set  B  that  attracts  points.  For  Equation 

(2.1) ,  a  set  B  attracts  bounded  sets  if  B  attracts  every 

bounded  subset  H  of  C.  A  set  B  is  stable  if,  for  any  neighbor¬ 

hood  U  of  B,  there  is  a  neighborhood  V  of  B  such  that 

T  (t)V  c  II  for  t  >  0.  A  set  B  is  asymptotically  stable  if  it  is 

stable  and  attracts  a  neighborhood  of  B. 

Definition  2.2.  The  positive  orbit  Y+(<{>)  through  <J>  is  the  set 
(T  _(t)$,  t  >  0).  A  function  h:  (-°°,0]  -*■  IR  is  said  to  be  a 

d  9  5 

backward  extension  of  4>  €  C  if  h„  =  $  and  T  (tlh.  =  h  .. 
-  u  a>g  i  t+Y 
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for  t  €  (-°°,0),  0  <  t  <  -t.  If  h  is  a  backward  extension  of  <J> , 


we  may 

define 

T  (t)4» 

for  t  5 

0  as 

Ta 

„(t)^ 

o 

V  I 

•M 

•C 

II 

A  ne 

gat i ve 

a»g 

a 

>g 

orbit 

Y"(<fr) 

is  the 

set  (To 

a»g 

(t)$, 

t 

<  0} 

defined  by 

some 

back- 

ward  extension 

of  4> . 

A  set  B 

in 

c 

is  invariant  if 

for 

every 

<t>  e  b, 

we  can 

define 

for 

t 

O 

V  1 

and  T  (t)<P  € 

a.g 

B  for 

t  6  ]R.  An  equilibrium  point  of  (2.1)  is  an  invariant  set  consisting 
of  a  single  point  <f> .  This  always  implies  is  a  constant 

function. 


A  periodic  orbit  of  (2.1)  is  an  invariant  set  which  is  a  closed 
curve;  that  is,  a  periodic  orbit  is  given  by  {p  ,  t  €  HR}  where 
p(t)  is  a  periodic  solution  of  (2.1). 

Let 

(2.2)  A  =  (4>  £  C:  T  (t)<|>  is  defined  and  bounded  for  t  <  0} 

a  >  8  a  >  g 

The  following  result  gives  some  of  the  fundamental  properties 

of  the  semigroup  T  (t)  and  the  set  A„  „  and  is  essentially 

a.g  a»g 

contained  in  Hale  [  3 ] ,  [  4 ] . 

Theorem  2.3.  If  Equation  (2.1)  is  point  dissipative,  then  Ag 
in  (2.2)  is  the  maximal  compact  invariant  set  of  Equation  (2.1)  and 

A„  _  is  uniformly  asymptotically  stable.  If,  in  addition, 

a»g 

Ta  (t)  is  one-to-one  on  Ag  then  Tg  ^(t)  is  a  continuous 

croup  on  A  .  Also,  if  T  (t)  takes  bounded  sets  into  bounded 

sets,  then  A  attracts  bounded  sets  of  C. 

- » -  a ,  g  - 

In  the  following  discussion,  S&  *  &  is  a  topological  space 

k 

which  is  also  a  subset  of  C  ( [  0 , 1  ]  JR)  *  C  (1RJR),  k  >  1,  with 
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Can,gn}  (a*g)  €  *  &  implying  that  an  -*■  a  uniformly  bn 

[0,1]  and  drgn/dxr  ->  drg/dxr,  0  <  r  <  k,  uniformly  on  compact  sets. 

Definition  2.4.  Stf  *  ^  is  uniformly  point  dissipative  if  there 

is  a  bounded  set  B  <=.  C  such  that  B  attracts  points  of  T  (t) 

a,g 

for  every  (a,g)  €  srf  *  %  . 

The  following  result  is  due  to  Cooperman  [  1]  (see,  also, 

Hale  [  4  ]) . 


Theorem  2.5.  If  x  5?  is  uniformly  point  dissipative,  then  A 

®  >  o 

is  upper  semicontinuous  in  ( a , g ) ;  that  is,  for  any  neighborhood  U 
of  A_  there  is  a  neighborhood  V  of  (a,g)  in  srf  *  %  such 

a  9  8 

that  A^  ^  c  U  for  (b,h)  €  V. 


Definition  2.6.  For  (a,g),  (b,h)  €  s/  *  &  ,  we  say 

equivalent  to  (b,h),  (a,g)  ~  (b,h),  if  there  is  a 

T:  A  A.  ,  such  that  XT  (t)  =  T,  ,  (t)x  on 

a,g  b,h  a,gv  1  b,!^ 

(a,g)  €  srf  *  &  is  structurally  stable  in  s/  *  9 
a  neighborhood  V  of  (a,g)  in  jrf  *  &  such  that 
for  every  (b,h)  €  V. 


(a ,g)  is 

homeomorphism 

A  .  An 
a  >g 

if  there  is 
(a ,g)  ~  (b,h) 


Definition  2.7.  An  equilibrium  point  a  of  (2.1)  is  hyperbolic  if 
the  solutions  of  the  characteristic  equation 

1 0  in 

(2.3)  X  +  g '  (ot )  a(-9)eAOd0  =  0 

}  -1 

have  negative  real  parts.  A  periodic  orbit  {pt,  t  G  1R}  is 
hyperbolic  if  the  linear  variational  equation 
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f° 

(2.4)  y(t)  +  a(-0)fi'(p(t+e))y(t+0)d0  =  o 

J  -1 

has  one  as  a  simple  characteristic  multiplier  and  no  other  character¬ 
istic  multiplier  on  the  unit  circle. 

If  p  has  period  t,  a  point  P  =  exp  At  is  a  characteristic 
multiplier  of  (2.4)  if  there  is  a  solution  of  (2.4)  of  the  form 
e^tq(t)  where  q  is  t-periodic  and  not  identically  zero  (see 
Hale  [  3  ]) . 

If  a  satisfies 


r° 

(2.5)  a(-0)d0  >  0 

1  -1 

then  an  equilibrium  point  a  is  hyperbolic  if  and  only  if  g'(«)  t  0, 
asymptotically  stable  if  g'(a)  <  0  and  unstable  if  g'(“)  >  0. 
Therefore,  the  set  of  all  g  for  which  the  equilibrium  points  are 
hyperbolic  is  a  residual  set  in  CrQR,lR)  for  each  fixed  a  satis¬ 
fying  (2.5). 

We  will  see  below  that  there  are  functions  a  for  which  there 
is  no  residual  set  of  g  for  which  the  periodic  orbits  are  hyperbolic. 
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3.  Convex  kernel.  In  this  section,  we  suppose  that  a  €  C  [-1,0], 


(3.1)  a(l)  =  0,  a(s)  >  0,  a(s)  <  0,  a(s)  >  0,  s  £  [0,1], 


(3.2)  G(x)  =  f  g  -*■  00  as  |  x |  -*• 

■>0 


The  following  result  may  be  found  in  Hale  [3,  p.  122] 


Theorem  3.1.  If  the  conditions  (3.1),  (3.2)  are  satisfied  and  the 
zeros  of  g  are  isolated,  then  every  solution  of  (2.1)  is  bounded, 
and 

( i )  If  there  is  an  s  such  that  a ( s )  >  0 ,  then  the  ^-limit 
set  of  any  orbit  of  (2.1)  is  an  equilibrium  point  of  (2.1);  that  is, 
a  constant  function  whose  value  is  a  zero  of  g .  An  equilibrium 
point  «  is  hyperbolic  if  and  only  if  g'(a)  f  0,  asymptotically 
stable  if  g'(a)  >  0  and  an  unstable  saddle  point,  if  g'(a)  <  0. 

(ii)  a(s)  =  0  (that  is,  a  is  linear)  then,  for  any  <J>  £ 

there  is  either  an  equilibrium  point  or  a  one-periodic  solution 
p  =  p(4>)  of  the  ordinary  differential  equation 

y  +  a(0)g(y)  =  0 

such  that  the  m-limit  set  of  the  orbit  through  <p  _is  (pt,  t  6  1R) 
For  any  zero  of  g,  let 


Ws(a)  =  { 4>  £  C:  T  (t)4>  -*■  the  constant  function  a  as  t  -*■  °°} 
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We  may  now  make  the  following  conjecture 

Conjecture  3.2.  I f  (3.1),  (3.2)  are  satisfied  and,  in  addition, 

the  zeros  ,  j  =  1,2,...,  o_f  g  are  simple  and  there  is  6  >  0 
such  that 


(3.3)  a (s )  >  0  for  s  €  (1-6,1) 

then  the  set 


W 


(Ws(a.): 


a  . 
J 


is  asymptotically  stable} 


is  dense  in  C . 

One  could  try  to  prove  this  using  an  idea  of  Henry  [  5  ] .  Let 

a  be  a  saddle  point  and  for  any  r  >  0 ,  t  >  0 ,  let  (a)  = 

(<J>  €  C:  T  (x )<J>  e  Ws(a),  |T  (t)$  -  a|  <  r).  The  set  w!(a)  is 

closed.  If  it  contains  a  ball  B  =  {4>^+ph:  |h|  =  1,  0  <  p  <  rQ>  , 

then  T  (t)B  has  been  flattened  so  that  it  is  contained  in  a 
3 » g 

submanifold  of  codimension  one.  In  particular,  9T  „  (t)  (<P  ,  +ph)/9p 

3  >  g  a 

c 

at  p  =  0  is  in  the  tangent  space  of  W  (a).  This  function 
Dx  T  D(t)h  is  the  solution  of  the  linear  variational  equation 

y  I  ’  “ 

y (t)  -  -  j  ^(-ejg'dg  gU^jMt) 

with  y^  =  h.  If  the  solution  operator  of  the  adjoint  of  this 
equation  is  one-to-one,  then  the  assertion  that  WT  r(a)  contains 
a  ball  is  false.  In  fact,  take  a  vector  n  /  0  such  that  <n,4>> 
for  all  <t>  in  the  tangent  space  to  Ws(cx)  at  Ta  g  ^ ^  j  . 
Integrating  the  adjoint  equation  for  a  function  z t  on  [0,t] 
with  zT  =  n,  one  obtains,  for  all  h, 
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(zQ,h)  =  (zT,h)  =  (n,D^Tajg(x)h)  =  0. 

Thus,  zQ  =  0.  Since  the  solution  operator  for  the  adjoint  equation 
is  one-to-one,  it  follows  that  n  =  0.  Consequently,  the  set 
Wsr  (a)  contains  no  ball,  which  implies  that  W®  is  closed  and 
nowhere  dense.  Thus  lVs(a)  has  no  interior.  This  will  prove  the 
conjecture . 

To  show  the  solution  operator  for  the  adjoint  equation  is  one- 

to-one,  it  is  sufficient  to  show  that  the  solution  operator  of  the 

linear  variational  equation  is  one-to-one  and  the  functions  which 

have  backward  extensions  are  dense  in  C. 

If  the  zeros  g  are  isolated  and  there  is  an  s  in  [0,1] 

such  that  a(s)  >  0,  then  Theorem  3.1  implies  every  solution  of 

(2.1)  approaches  an  equilibrium  point  of  g.  If  we  assume,  in 

addition,  that  there  are  only  a  finite  number  of  zeros  of  g,  then 

the  method  of  proof  of  Theorem  3.1  in  Hale  [  3]  shows  that  A  is 

a  9  8 

compact  and  attracts  bounded  sets  of  C.  Also, 


Aa,g  ' 


where  is  a  zero  of  g  and  Wu(Oj)  is  the  unstable  manifold 

of  dj .  The  unstable  manifold  of  any  zero  has  dimension  <1  as 

is  easily  seen  from  Equation  (2.5).  Thus,  is  the  union  of  a 

finite  number  of  manifolds  of  dimension  <  1.  If  each  ct is 

hyperbolic  and  the  w-limit  set  of  each  unstable  manifold  Wu(a^) 

of  an  unstable  equilibrium  point  is  a  stable  equilibrium  point, 

then  A  has  the  one-dimensional  structure  shown  in  Figure  1. 
a ,  g 
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• - — *-5» - • - «Sr- 


Figure  1 . 

Suppose  srf  *  &  is  a  uniformly  bounded  dissipative  set  containing 

(a,g).  To  see  how  to  construct  such  a  set,  consult  the  proof  of 

Theorem  3.1  in  Hale  [3].  We  can  show  that  (a,g)  is  structurally 

stable  in  srf  *  &  .  In  fact,  Theorem  2.5  implies  ^  if  upper 

semicontinuous  at  (a,g).  Also,  there  is  a  neighborhood  U  of 

A„  such  that  the  flow  is  given  as  in  Figure  2- 
a,g  6  b 


Figure  2. 


This  neighborhood  U  can  be  chosen  to  be  invariant  for  all  (b,h)  £  V. 

Thus  the  new  A,  , ,  (b,h)  £  V  must  have  the  same  structure  as  A 

o  ,  n  a  ,  g 

We  state  this  as 
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Theorem  3.3.  Suppose  there  is  an  s  in  (0,1)  such  that  a(s)  >  0 
and  there  are  only  a  finite  number  of  equilibrium  points  of  (2.1), 
each  of  which  is  hyperbolic.  If  the  u>-limit  set  of  an  unstable 
manifold  is  a  stable  equilibrium  and  *  &  is  a  uniformly  bounded 
dissipative  set  containing  (a,g),  then  (a,g)  is  structurally 
stable  in 

One  can  now  ask  the  following  interesting  question. 

Question  3.4.  Suppose  conditions  (3.1),  (3.2)  are  satisfied  and 
a  is  a  fixed  function  satisfying  a(s)  >  0  for  s  G  (0,1).  Let 
-^k  =  which  have  exactly  2k+l  zeros  which  are  all  simple). 

How  many  different  equivalence  classes  are  in 

If  k  =  0,  that  is,  g  has  only  one  zero  «,  then  A  =  {« } , 

the  constant  function  a  and  all  g  G  ^  are  equivalent. 

If  k  =  1,  then  g  G  ^  has  three  simple  zeros,  <  “2  <  a3» 

with  >a3  asymptotically  stable  and  a ^  a  saddle  point.  Since 

the  unstable  manifold  at  a~  is  smooth,  one  dimensional  and  A 

i-  a  >g 

is  uniformly  asymptotically  stable,  it  follows  that  A  is  a  one 

a ,  g 

dimensional  manifold  with  boundary  points  ai»a3*  Again,  all  ele¬ 
ments  of  are  equivalent. 

The  topological  structure  of  A  is  not  understood  for  the 

case  when  g  has  given  zeros  ai  <  a2  <  °3  c  a4  <  a5  •  If  we  consider 

g  depending  on  some  parameters  M  with  g(x,W)  =  x^  at  W  =  0, 

then  A  (y)  is  {0}  for  y  =  0.  There  is  a  center  marifold  in 
a  t  g 

a  neighborhood  of  x  =  0  which  is  one  dimensional  and  smooth  in  y. 
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As  vi  varies  near  u  =  0,  g(x,p)  can  have  as  many  as  5  zeros 

which  must  lie  on  this  center  manifold.  Thus,  A  (u)  is  either 

a ,  g 

a  point  or  a  one-dimensional  manifold  with  boundary  for  M  small. 
As  p  increases,  it  is  conceivable  that  the  topological  structure 
of  Aa  changes.  Let  us  give  some  intuitive  reasons  for  why  this 
is  possible.  The  author  wishes  to  acknowledge  conversations  with, 
John  Mallet-Paret  which  were  of  great  assitance  in  the  remaining 
discussion  of  this  section. 

Suppose  g  has  five  zeros  and  the  general  shape  shown  in 
Figure  3.  If  a  is  strictly  convex  and  very 


Figure  3.  Figure  5  . 


close  to  the  (5-function  at  zero,  then  A  is  shown  in  Figure  4. 

a  >  8 

The  complete  flow  near  A  is  shown  in  Figure  5  with  the 

a,g 

exponential  decay  toward  Afl  being  very  rapid  and  much  greater 

than  the  convergence  of  the  flow  on  A  toward  6. 

a,g 

If  g'({3)  =  a>  then  the  flow  near  (3  is  determined  by  the 
roots  of  the  characteristic  equation 
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(3.4) 


(° 

■“  J-i 


a(-6)eX0d9 


One  can  show  there  is  an  and  a  strictly  convex  function  a^ 

with  the  property  that  the  roots  of  this  equation  with  maximum  real 
part  corresponds  to  a  double  root  root.  Furthermore,  in  any  neighbor¬ 
hood  of  (aQ,ap)  there  is  an  and  a  strictly  convex  function  a 

such  that  Equation  (3.4)  has  a  pair  of  complex  conjugate  roots 
X0,X0  with  maximum  real  parts.  For  this  a  and  the  corresponding 
g  and  from  the  fact  that  the  flow  near  Ag  is  shown  in  Figure  5 
with  the  stable  manifolds  transversal  to  A  of  codimension  one, 

a 

the  flow  near  Ag  is  governed  by  what  happens  on  a  two-dimensional 
manifold  Mg  determined  by  the  invariant  manifolds  corresponding 

to  the  roots  This  manifold  is  constructed  by  finding  a  local 

invariant  manifold  near  0  which  is  tangent  to  the  invariant  sub¬ 
space  of  the  linear  approximation  near  0  generated  by  the  eigen¬ 
values  A  ,X_.  Then  find  M  by  backward  extension  along  the 
0  0  a  ,g 

flow.  On  M  the  flow  is  shown  in  Figure  6. 
a,g 


Figure  6. 


Figure  7. 
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If  it  is  possible  to  show  that  one  can  break  the  connection 
between  Y,6  by  a  small  variation  in  a,  then  the  new  picture  will 
be  as  in  Figure  7.  This  is  the  point  that  the  author  has  been 
unable  to  verify. 

Assuming  that  the  previous  step  can  be  done,  we  have  changed 
the  orientation  of  the  points  on  ;  namely  from  Figure  8a) 

to  Figure  8b).  Thus,  they  are  in  different  equivalence  classes. 


Figure  8a) 


Figure  8b) 


A  small  change  in  g  to  make  the  double  zero  either  disappear 
or  give  two  simple  ones  gives  the  distinct  structurally  stable  systems 
shown  in  Figure  9.  We  draw  the  spirals  as  arcs  since  they  are  the 


Figure  9a) 


Figure  9b) 


Figure  9c) 
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same  under  homeomorphism. 

Now  let  us  vary  g  some  more  to  attempt  to  go  from  the 
structurally  stable  system  in  Figure  9c)  to  the  one  in  Figure  9b) 
and  keep  the  zeros  of  g  simple.  The  only  possible  way  to  do  this 
is  to  have  a  saddle  connection;  namely,  have  a  function  g  for 
which  the  saddle  6^  is  connected  to  Y  with  the  configuration 
shown  in  Figure  10. 


Figure  10. 

Conjecture  3.5.  There  is  a  saddle  connection  for  some  (a,g) 
with  a  strictly  convex  and  g  having  five  simple  zeros.  The  set 
of  (a,g)  for  which  this  is  true  is  not  generic. 

Theorem  3.1  can  be  generalized  to  an  n-vector  system  of  the 

form 

[0 

(3.5)  x ( t )  =  -  a(-6)grad  G(x(t+0))d0 

1  -1 

where  x  e  lRn,  G:  lRn  -*■  1R.  The  proof  will  be  the  same  as  in  Hale  [  3  ]. 
A  proof  or  counterexample  to  the  following  conjecture  would  be 


interesting: 


Conjecture  3.6.  I f  a(l)  =  0,  a(s)  >  0  for  0  <  s  <  1,  then  the 
set  of  G  for  which  Equation  (3.5)  is  structurally  stable  is 
residual  in  the  set  of  C ^  functions  from  lRn  _to  1R  with  the 
Whitney  topology. 

Let  us  now  consider  in  more  detail  the  case  when  the  kernel 
in  Equation  (2.1)  is  a  linear  function.  Consider  the  equation 

fO 

(3.6)  x(t)  =  -\i  j  ao(-9)g(x(t  +  0))d8,  aQ(s)  =  (1-s),  U  >  0 

and  assume  that  the  zeros  of  g  are  isolated.  From  part  (ii)  of 
Theorem  3.1,  we  know  that  the  w-limit  set  of  every  solution  of  (3.6) 
is  either  an  equilibrium  point  of  a  periodic  orbit  generated  by  a 
1-periodic  solution  of  the  ordinary  differential  equation 

(.3.7)  y  +  ug (y )  *  0.  . 

As  we  shall  see  below,  the  structure  of  A  can  be  very 

a  ,g 

complicated  and  is  not  very  well  understood.  However,  the  following 
interesting  fact  is  true. 

Proposition  3.6.  For  Equation  (3.6),  no  periodic  orbit  can  be 
uniformly  asymptotically  stable  relative  to  perturbations  in 
A.._  _.  If  A.._  „  contains  a  periodic  orbit,  then  (pan,g) 

H«Qf8  »6  " 

cannot  be  structurally  stable. 

wa0,g 


Proof:  Suppose  A, 


contains  a  periodic  orbit  y  which  is 
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uniformly  asymptotically  stable.  Then,  for  any  neighborhood  U  of 
Y  there  is  a  neighborhood  V  of  Y  and  a  neighborhood  W  of 
Map  in  the  C°-topology  such  that,  for  any  a  €  W,  <f>  €  V,  the 
positive  orbit  Y+(<i>)  of  (2.1)  belongs  to  U.  In  particular,  one 
can  choose  the  neighborhood  U  of  Y  and  W  of  paQ  so  that  U 
contains  no  equilibrium  point  of  (2.1)  for  any  a  e  W. 

In  the  neighborhood  IV  of  ua^  there  exists  a  strictly 
convex  function  a.  Thus,  A  contains  no  periodic  orbits  and 

a,g  ,  p 

every  solution  of  (2.1)  approaches  an  equilibrium  point.  This 

contradicts  the  fact  that  some  positive  orbits  remain  in  U.  This 

proves  no  periodic  orbit  can  be  uniformly  asymptotically  stable. 

The  fact  that  (an,g)  cannot  be  structurally  stable  if  A 

u  ya0’g 
contains  a  periodic  orbit  follows  by  choosing  a  strictly  convex 

function  near  uaQ.  This  proves  the  proposition. 

To  understand  more  about  the  possible  topological  structures 

for  A^a  g,  we  consider  some  special  cases.  Suppose 


yg(y)  >  0,  y  f  0,  g'(0)  =  1, 


If  p(t,G)  is  the  solution  of  the  equation 


g(y) 


with  p(0,b)  =  b  >  0,  p ( 0 , b )  =  0,  then  p(t,b)  has  period 

w(b)  2 it  as  b  -*•  0 .  We  assume  that  u> ( b )  is  a  strictly  monotone 

function  of  G  for  b  >  0. 

1/2 

The  solutions  of  (3.7)  are  then  given  by  p(u  '  t,b)  and  have 

1/2 


period  w(b)/y 
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The  linear  variational  equation  around  x  =  0  is 


(3.8) 


■  -  l-i 


x ( t )  =  -v  (l+0)x(t+0)d0 


for  which  the  characteristic  equation  is 


(3.9) 


(° 

"  1-1 


(1+9  )e 


X0 


or 

(3.10)  X  +  U(A-l+e'X)A'2  =  0. 

Equation  (3.7)  is  a  special  case  of  Equation  (2.1)  and  so  Theorem  3.1 
may  be  applied.  Thus,  we  see  that  every  solution  of  (3.7)  approaches 
zero  as  t  -*■  °°  if 


(3.11)  M  *  4Tr2k2,  k  =  1,2,3 . 

This  implies  the  roots  of  Equation  (3.9)  have  negative  real  parts 
if  y  satisfies  (3 . 10) . Consequently ,  the  solution  x  =  0  of 
Equation  (3.6)  is  uniformly  asymptotically  stable  if  y  satisfies 
Relation  (3.10). 

2  2 

For  u  =  4tt  k  ,  Equation  (3.9)  has  the  purely  imaginary  solutions 

+  A,  \  =  2 Ttk i  with  all  other  solutions  having  negative  real 

2  2 

parts.  Also,  the  solution  A  (y )  with  \(4it  k  )  =  2irki  satisfies 
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dX  =  i 
Hy  ”  &tik 


24TT4k 


-) 


at  y  =  4^  2k2 
at  y  =  4^2k2. 


Thus,  Re  X  behaves  as  a  quadratic  function  and  Im  X  as  a  linear 

2  2 

function  in  a  neighborhood  of  y  =  4^  k  . 

We  can  begin  to  use  this  information  to  analyze  the  set  A 

yaQ,g 

2 

for  Equation  (3.6),  If  0  <  u  <  4^  ,  there  are  no  1-periodic  solutions 

of  Equation  (3.7)  and  every  solution  approaches  zero  as  t  -*■  ®. 

2  2 

Thus,  A^a  =  (0).  If  4^  <  y  <  16^  ,  there  is  exactly  one 

1-periodic  orbit  Y^  in  A^g  generated  by  the  1-periodic  solution 

of  Equation  (3.7).  Also,  {0}  c  A  and  the  solution  x  =  0  is  uni- 

uaQ,g 

formly  asymptotically  stable  by  the  previous  analysis  of  Equation  (3.7). 

For  y  very  close  to  zero,  there  is  a  two  dimensional  center  manifold 

M'  which  contains  A,,  .  Since  A,,  is  uniformly  asymptotically 

^  Ma0,g  ^a0’g 
stable,  it  follows  that  A  is  the  two  manifold  consisting  of  the 

VJa0,g 

orbit  Y '  and  its  interior  relative  to  M' .  The  flow  on  A  is 

y  y  ua0,g 

given  as  in  Figure  11.  The  orbit  yJ  in  the  center  manifold 


Figure  11. 


Figure  12. 
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is  asymptot i cal ly  stable  from  the  outside  and  unstable  from 

the  inside.  The  orbit  Yj  has  an  “-dimensional  stable  manifold  which 

is  transversal  to  so  that  locally  the  flow  is  similar  to  the  one 

shown  in  Figure  12.  Everything  inside  the  ’’cylinder"  through  y^ 

consisting  of  the  stable  manifold  transverse  to  approaches  zero 

and  everything  outside  approaches  Y^ .  If  we  vary  y  in  the  interval 

4*  <  y  <  16*  ,  the  same  picture  will  prevail  since  0  and  Y^  must 

be  the  attractor  for  all  orbits.  Thus,  A  is  a  "disk"  for 

ua0,g 

4*2  <  y  <  16*2  and  is  given  by  A  =  {0}  U  WU(Y'). 

2  ya0’g 

At  y  =  16*  ,  there  is  a  solution  of  (3.9)  on  the  imaginary  axis 

2  2 

given  by  X  =  4tt i  and  for  16tt  <  y  <  36tt  ,  there  is  a  unique 

periodic  orbit  yJ  of  (3.6)  of  period  1/2  (and  therefore  of  period  1) 

2 

generated  by  the  (1/2 -periodic  solution  of  (3.7).  For  16n  <  y  < 

36*  ,  the  origin  is  uniformly  asymptotically  stable  from  the  above 

analysis  of  Equation  (3.8).  Again,  there  is  a  two-dimensional  center 
2 

manifold  through  x  =  0  which  is  transversal  to 

{0}  U  WU(Y'  ?)  and  asymptotically  stable.  The  orbit  Y. 

16*^  M 

2  2 
must  be  on  .  Reasoning  as  before,  we  have  for  y  >  16*  but 

close  to  16*2  that  A  n  M2  =  (0}uWu(Y,f).  Since  Wu(y2)  is 

y“g»gy  y  y 

transversal  to  Wu(y'),  the  set  A  is  two-dimensional  and  is 

y'  yaQ’S 

the  union  of  a  finite  number  of  manifolds.  Locally  near  0,  the 

set  A  is  like  two  disks  in  4-dimensional  space  intersecting 

yaQ  >  g 

at  zero. 

2 

As  y  passes  through  36*  ,  another  priodic  orbit  of  period 
1/3  appears  and  something  similar  happens,  but  the  geometry  becomes 
more  complicated. 
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If  the  function  g  in  (3.6)  has  three  simple  zeros,  say 

a  <  6  <  y,  then  3  must  always  be  an  unstable  saddle  point  with 

Wu(8)  one-dimensional.  The  equilibrium  points  a,Y  will  be  uni- 

2  2 

formly  asymptotically  stable  as  long  as  Mg' (a)  t  4^  k  , 

2  2 

U|>'  (8)  /  4"  k  for  all  k  =  1,2,...  .  The  phase  portrait  for  the 
second  order  equation  (3.7)  is  shown  in  Figure  13.  If  there  is  no 


Figure  14b) 


periodic  orbit  of  period  1,  then  A,,„  is  a  one-dimensional  mani- 

.  ^a0,g 

fold  with  boundary  points  being  the  constant  functions  <*,Y. 

Suppose  now  that  there  is  no  periodic  solution  of  period  1  whose 

orbit  encircles  all  three  equilibrium  points  and  only  one  encircling 

each  of  the  points  (a»0)»(Y,0).  Arguing  as  for  the  previous  case, 

one  sees  that  A  has  the  structure  shown  in  Figure  14a)  with 

Wa0>8 

the  flow  on  A  shown  if  Figure  14b).  Each  of  the  periodic 

yaQ,g 

orbits  is  unstable  from  the  inside  on  A^  ,  Suppose  now  that 

nag ,  g 
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there  is  no  1-periodic  orbit  of  (3.7)  encircling  one  equilibrium 

point  and  exactly  one  1-periodic  orbit  y  encircling  all  three 

equilibriums.  We  do  not  know  the  structure  of  A,  but  con- 

Pa0>R 

jecture  that  it  is  a  two-dimensional  manifold  with  boundary  Y, 

the  equilibrium  points  ct,Y  are  asymptotically  stable  and  3,Y 

are  unstable  on  A, 

Ma0,g 

By  choosing  different  functions  g,  one  can  certainly  make  the 


set  A  have  a  large  number  of  periodic  orbits  and  equilibrium 

^aQ  *  8 

points.  The  way  these  periodic  orbits  and  equilibrium  points  fit 


together  in  the  flow  on  A 


Va0,g 


is  not  known.  However,  from 


Proposition  3.5,  we  know  that  no  periodic  orbit  in  A  f  can  be 

a0’r 

uniformly  asymptotically  stable. 

Fron  the  above  analysis,  we  see  that  Equation  (2.1)  with  a(s) 
strictly  convex  has  a  very  nice  structure  which  appears  to  be 
structurally  stable  for  a  residual  set  of  functions  g.  Also,  one  can 
show  that  (a,g)  structurally  stable  for  -  a  strictly  convex  will 
imply  that  (a,g)  is  structurally  stable  if  g  is  close  to  g  and 
a  is  close  to  a;  that  is,  the  kernel  also  can  be  cahnged  slightly. 

If  aQ(s)  is  linear,  it  is  probably  true  that  (aQ,g)  is 
structurally  stable  for  a  residual  set  of  g.  However,  (a0,gQ) 
structurally  stable  subject  to  variations  in  g  does  not  imply  that 
(a,g)  will  be  structurally  stable  if  one  changes  a^  slightly  to  a 
new  function  a  (Proposition  3.5).  In  fact,  the  structure  of  A 

a  >  g 

can  be  drastically  different  for  a  strictly  convex  function  a  which 

2 

is  close  to  Sq  in  the  C  -topology.  In  fact,  all  periodic  orbits 
in  A,.  will  disappear.  Thus,  it  becomes  interesting  to  analyze 

the  behavior  of  the  solutions  of  (2.1)  for  an  arbitrary  function  a 
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2 

close  to  the  linear  function  Pa()  in  the  C  -topology.  This  is 
the  topic  disuccsed  in  the  next  section. 
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4,  Bifurcation  near  linear  kernels.  In  this  section,  we  consider 
equation  (2.1)  with 


(4.1) 


xg(x)  >  0  for  x  0,  g'(0)  =  1, 


i(.x)  = 

i 


as  I  x  I  -*■«>, 


where  g  is  continuous  together  with  all  derivative  up  through 

0  2 

order  five  and  a  is  in  a  C  -neighborhood  of  4 rr  aQ,  aQ(s)  =  1  -  s 
We  also  suppose  the  period  oo(b)  of  the  solution  through  (b,0), 
b  >  0,  of  the  equation 


(4.2) 


satisfies 


x  +  g(x)  =  0 


(4.3) 


<•>’  (b)  f  0. 


Under  an  additional  conditional  on  g  which  will  be  specified 

later,  we  determine  the  behavior  of  the  solutions  of  (2.1)  in  a 

neighborhood  of  x  =  0  for  all  a  in  a  neighborhood  of  the  bi- 

2 

furcation  point  4*  a^.  The  analysis  is  based  on  the  application 
of  the  method  of  Liapunov-Schmidt  for  periodic  orbits  near  x  =  0 
together  with  properties  of  the  bifurcation  function  from 
deOliveira  and  Hale  [2]. 

We  need  the  following  lemma. 


2  0 

Lemma  4.1.  There  is  a  neighborhood  U  of  4^  aQ  in  the  C  -topology, 
a  6  >  0  and  an  analytic  function  :  U  ■+  C  such  that  X*(4^an)  =  2*i 
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and,  for  every  a  £  U,  the  equation 

fO  \  A 

(4.4)  X  +  a(0)e  d6  =  0 

J  -1 

has  exactly  one  solution  in  each  of  the  circles  |X  ±  2^i\  <  6 
given  respectively  by  X*(a),^  (a)  and  all  other  solutions  with  real 

parts  <  -6.  Furthermore,  if 

r"  -  {a  €  U:  Re  X*(a)  <  0} 

T°  =  {a  €  U:  Re  X*(a)  =  0} 

T1  =  (a  £  U:  Re  X*(a)  >  0} 

then  each  of  these  sets  is  nonempty  and  r°  is  a  submanifold  of 

codimension  one. 


Proof.  I f 

{0  \  a 

F (X , a)  =  X  -  a(8)e  d0 

1  -1 

then  F(2Tii  ,4Tr2a0)  =  0  and  3f/9X  =  1  -  8"2  at  (2ni  ,4n2a0) .  The 

Implicit  Functfan  Theorem  implies  the  existence  of  a  function 

X*(a)  e  analytic  in  a  neighborhood  of  aQ  with  X*(aQ)  =  2 tt i . 

The  other  properties  of  X*  follow  essentially  fromRouche's  Theorem. 

To  show  r°  has  codimension  one,  consider  the  family  of 
2 

functions  4n  (aQ  +  vb0),  b ^ ( s )  =  s(l-s),  v  e  IR.  Then  the  derivative 
of  x*  ( 4 tt 2  (a^*  vbg)  with  respect  to  v  at  v  =  0  is  easily  seen 


to  satisfy 


SHv-o 


- 4tt2  j°  o  Cl+0)e27ri9d0 


Thus , 


9  Re  X' 


2  f° 

=  -47T  I 

>=0  M 


9  (l-9)cos  2iT9d9  >  0. 


This  shows  that  r°  has  codimension  one  and  also  that  T  ,r  +  arc 
not  empty.  This  proves  the  lemma. 


Remark  4.2.  We  knew  from  the  previous  section  that  Re  X*(a)  <  0 
if  a  €  U  is  strictly  convex.  The  above  proof  shows  that 
Re  X*(a)  >0  if  a  G  U  is  of  the  form  a  =  4rcZaQ  +  b  where  b 


is  strictly  concave,  b(0)  =  b ( 1 )  =  0. 


For  a  =  4tt  aQ,  the  characteristic  equation  for  the  linear 

variational  equation  around  zero  has  two  purely  imaginary  roots. 

2 

For  a  near  4n  a q,  and  a  neighborhood  W  of  zero  let  B(r,a,g) 

be  the  scalar  bifurcation  function  obtained  by  applying  the  usual 

method  of  Liapunov-Schmidt  for  the  periodic  solutions  of  (2.1)  in 

2 

W  which  for  a  =  4u  a0  are  equal  to  r  cos  2* t  (see,  for 
example,  dcOliveira  and  Hale  [  2]).  This  function  has  the  property 
that  the  periodic  solutions  of  the  type  specified  are  in  one  to  one 
correspondence  with  the  nonnegative  zeros  of  B(r,a,g).  Further¬ 
more,  the  stability  properties  of  the  periodic  solution  correspond¬ 
ing  to  a  zero  Tq  of  B(r,a,g)  when  restricted  to  a  center 
manifold  at  x  =  0  are  the  same  as  the  stability  properties  of  the 
equilibrium  point  rn  of  the  scalar  equation 


(4.5) 


r  =  B (r ,a , g) 


(?ee  deOliveira  and  Hale  (2]).  The  function  B(r,a,g)  is  an  odd 
function  of  r  and  has  five  continuous  derivatives.  Let 

(4.6)  B (r , a ,g)  =  «1(a,g)r  +  “3(a,g)r3  +  a5(a,g)r5  +  o(|r|5) 

as  r  ->  0. 

If  X*(a)  is  the  function  given  in  Lemma  4.1,  the  manner  in 
which  the  bifurcation  function  is  constructed  implies  that 

ai(a,g)  =0  if  and  only  if  Re  X*(a)  =  0 

(4-7) 

sign  “1(a,g)  =  sign  Re  X  (a) 

Thus,  «1(4Ti2a0,g)  =  0. 

Let  a3  =  <*3  C 4 n" ^a^ , g) .  We  now  show  that  =  0.  Since  the 

2 

solution  x  =  0  of  Equation  (2.1)  for  a  =  4^  is  asymptotically 

2 

stable,  it  follows  that  the  zero  solution  of  (4.5)  for  a  =  4*  a0 
is  asymptotically  stable.  Thus,  °3  <  0.  If  «3  <  0,  then  the 

2 

function  B(r,a,g)  for  a  varying  in  a  small  neighborhood  of  4 n 
can  have  only  three  possible  shapes  in  a  neighborhood  of  r  =  0  as 
shown  in  Figure  15. 

Now  consider  the  special  function  a  =  paQ  for  p  varying  in 

2  2  2 

small  neighbhood  of  4 it  ,  |  ia - 4 tt l\  <  e.  If  p  /  4*  ,  we  know  that 

Re  X*(pap)  <  0  from  the  previous  section.  Thus,  “jtuag.g)  <  0 
2 

if  p  f  4it  .  For  definiteness,  suppose  w'(b)  <0  in  (4.3).  Then, 


Figure  15. 


the  analysis  in  the  previous  section  implies  there  is  an  e  >  0  and 
a  sufficiently  small  neighborhood  W  of  x  =  0  such  that  there  is 

7 

no  periodic  orbit  of  (2.1)  in  W  for  0  <  p  -  4 n  <  e  and  a  unique 

2 

periodic  solution  y^  in  W  for  -e  <  p  -  4n  <  0 ,  ■+  0  as 

2 

U  -*■  4^  and  y  is  unstable  in  A„  .  Thus  on  a  center  manifold, 
u  a,g 

Yp  is  stable  from  the  outside  and  unstable  from  the  inside.  This 
means  the  same  is  true  for  the  corresponding  equilibrium  point  r^ 
Equation  4.5.  Consequently,  the  bifurcation  function  must  have  the 
shape  shown  in  Figure  16.  This  is  impossible  if  <  o.  Thus, 


Figure  16. 
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Thus,  “j  =  0. 

Again,  the  stability  of  the  zero  solution  of  Equation  (2.1)  for 
2  2 

a  =  4 w  a^  implies  that  <*^.(471  ag,g)  <  0.  We  make  the  hypothesis 
that 

(4.8)  “5(g)  d=f  a5(4Tr2a0,g)  <  0. 

This  implies  that 

(4.9)  B(r,4*2a0,g)  =  «°(g)r5  +  o(r5),  <  0. 

We  have  not  made  the  computations  (which  would  be  extremely  complicated) 

to  obtain  the  constant  a[j(g).  However,  it  certainly  seems  plausible 

that  the  set  of  g  for  which  a^(g)  <  0  is  open  in  the  space 

C^(U^IR)  for  a  given  bounded  neighborhood  U  of  x  =  0. 

2 

If  B(r ,a,g)  =  rP(r  ,a,g),  then 

(4.10)  P(P,a,g)  =  0l1  (a,g)  +  a  3  ( a ,  g )  P  +  a5(a,g)P2  +  o(P2) 

as  p  -*■  0.  This  function  has  a  unique  maximum  n(a,g)  in  a 
neighborhood  U  of  a  =  41r2ag  which  occurs  at  a  value  p*(a,g)  and 
n(4n2a0,g)  =  0.  If 
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(4.11) 


SN°  =  {a  G  U:  n(a,g)  =  0,  P*(a,g)  >  0) 
SN+C_)  =  {a  €  U:  n(a,g)  >  (<)0,  P*(a,g)  >  0}. 


One  can  show  that  every  tangent  vector  to  SN^  is  a  tangent  vector 
to  r  .  We  suppose  that 

(4.12)  SN^  is  a  submanifold  of  codimension  1,  S+  f  <t> . 


It 

an  open 
We 


is  possible  to  show  that  hypothesis  (4.12)  is  satisfied 
set  of  g  €  CS(W,1R)  for  a  bounded  neighborhood  W  of 
can  now  prove  the  following  result. 


for 

zero 


Theorem  5.3.  With  hypotheses  (4.1),  (4.3),  (4.8)  and  (4.12),  there 
is  a  neighborhood  U  of  4^  aQ  in  the  C  -topology  and  a  neighbor¬ 
hood  W  of  x  =  0  such  that  U  is  subdivided  into  regions  as 
shown  in  Figure  17,  the  set  A„  _  is  a  disk  for  each  a  €  U  with 

- -  -  a ,  g  - 

boundary  being  a  periodic  orbit  and  the  flow  on  a  two  dimensional 

manifold  in  A_  _  is  shown  in  Figure  17. 

~  a , g 


II 
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Proof.  If  a  €  r  +  ,  then  01  -j  C.a , g )  >  0,“3Ca,g)  <  0  implies  that 

P(P»a>g)  in  (4.10)  has  a  unique  positive  zero.  Thus,  there  is  a 

unique  periodic  solution  in  a  small  neighborhood  of  zero  and  it  is 

asymptotically  stable  as  shown  in  the  flow  for  r+.  This  shows  that 

SN  c  r*.  By  hypothesis  (4.12)  and  the  fact  that  the  stability 

properties  of  the  priodic  orbits  are  determined  by  (4.5),  we  have  that 

the  flow  on  Sf/*  is  the  one  shown  in  Figure  17.  Also,  the  flow  in  the 

other  two  regions  must  be  one  of  those  shown  in  Figure  17.  We  only 

need  to  verify  that  the  regions  are  ordered  as  shown.  In  the  proof  of 

2 

Lemma  4.1  we  showed  that  the  curve  4it  (a0+vbQ) ,  bg(s)  =  s(l-s)  was 
transversal  to  at  v  *  0.  For  v  <  0  this  function  is  in  Y 

and  stictly  convex.  Thus,  the  origin  is  uniformly  asymptotically 
stable.  This  proves  that  the  flow  in  Region  II  in  Figure  17  is  the 
one  that  is  depicted.  This  proves  the  theorem. 
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